arXiv: 1506.01155vl [math.PR] 3Jun2015 


On the first-passage time of an integrated 
Gauss-Markov process 

Mario Abundo* 


Abstract 

It is considered the integrated process X{t) = x + ff^Y(s)ds, where y(t) is a Gauss- 
Markov process starting from y. The first-passage time (FPT) of X through a constant 
boundary and the first-exit time of X from an interval (a, b) are investigated, general¬ 
izing some results on FPT of integrated Brownian motion. An essential role is played 
by a useful representation of X, which allows to reduces the FPT of X to that of a 
time-changed Brownian motion. Some explicit examples are reported; when theoret¬ 
ical calculation is not available, the quantities of interest are estimated by numerical 
computation. 
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1 Introduction 

First-passage time (FPT) problems for integrated Markov processes arise both in theoretical 
and applied Probability. For instance, in certain stochastic models for the movement of a 
particle, its velocity, D(t), is modeled as Ornstein-Uhlenbeck (OU) process, which is indeed 
suitable to describe the velocity of a particle immersed in a fluid; as the friction parameter 
approaches zero, Y {t) becomes Brownian motion Bt (BM). More generally, the particle 
velocity Y{t) can be modeled by a diffusion. Thus, particle position turns out to be the 
integral of D(t), and any question about the time at which the particle first reaches a given 
place leads to the FPT of integrated Y[t). This kind of investigation is complicated by the 
fact that the integral of a Markov process such as Y{t), is no longer Markovian; however, the 

two-dimensional process ^(t) = y(s)(is, is Markovian, so the FPT of integrated 

Y(t) can be studied by using Kolmogorov’s equations approach. The first apparition in the 
literature of ^(t), with Y{t) = Bt , dates back to the beginning of the twentieth century 
(see [23]), in modeling a harmonic oscillator excited by a Gaussian white noise (see [2l] and 
references therein). 

The study of Y{s)ds has interesting applications in Biology, in the framework of diffu¬ 
sion models for neural activity; if one identifies Y[t) with the neuron voltage at time t, then 
j f*Y(s)ds represents the time average of the neural voltage in the interval [0,t]. Moreover, 
integrated Brownian motion arises naturally in stochastic models for particle sedimentation 
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in fluids (see |2I])- Another application can be found in Queueing Theory, if Y(t) represents 
the length of a queue at time t; then f*V(s)ds represents the cumulative waiting time ex¬ 
perienced by all the “users” till the time t. Furthermore, as an application in Economy, one 
can suppose that Y(t) represents the rate of change of a commodity’s price, i.e. the current 
inflation rate; hence, the price of the commodity at time t is X{t) = X(0) -|- Y(s)ds. 
Finally, integrated diffusions also play an important role in connection with the so-called 
realized stochastic volatility in Finance (see e.g. i. 1151. lEl). 

FPT problems of integrated BM (namely, when Y(t) = B^) through one or two bound¬ 
aries, attracted the interest of a lot of authors (see e.g. [TO], [I7j, [21], [25], [26], [28], [34] for 
single boundary, and m, 1311. 1321 for double boundary); the FPT of integrated Ornstein- 
Uhlenbeck process was studied in [10], [29]. Motivated by these works, our aim is to extend 
to integrated Gauss-Markov processes the literature’s results concerning FPT of integrated 
BM. 

Letm(t), h 2 (t) be (PAfunctions of t > 0, such that h 2 (t) 7 ^ 0 and p(t) = hi{t)/h 2 {t) 

is a non-negative and monotonically increasing function, with p( 0 ) = 0 . 

If B{t) = Bt denotes standard Brownian motion (BM), then 

Y{t) = m{t) + h 2 {t)B{p{t)), t>0, (1.1) 

is a continuous Gauss-Markov process with mean m{t) and covariance c{s,t) = hi{s)h 2 {t), 
for 0 < s < t. 

Throughout the paper, Y will denote a Gauss-Markov process of the form fll.ip . starting 
from y = m( 0 ). 

Besides BM, a noteworthy case of Gauss-Markov process is the Ornstein-Uhlenbeck (OU) 
process, and in fact any continuous Gauss-Markov process can be represented in terms of a 
OU process (see e.g. 1351). 

Given a continuous Gauss-Markov process Y, we consider its integrated process, starting 
from X( 0 ) : 

X{t) = A(0) + [ Y{s)ds. (1.2) 

Jo 

For a given boundary a, we study the FPT of X through a, with the conditions that X(0) = 
X < a and Y ( 0 ) = y, that is: 

Taix,y) = inf{t > 0 : X{t) = a|X(0) = x,Y{0) = y}; (1.3) 

moreover, for b > a and x G {a,b), we also study the first-exit time of X from the interval 
(a, b), with the conditions that 2 f( 0 ) = x and Y ( 0 ) = y, that is: 

Ta,b{x,y) = inf{t >0: X{t) ^ (a,5)|X(0) = x,U(0) = y}. (1.4) 

In our investigation, an essential role is played by the representation of X in terms of BM, 
which was previously obtained by us in [T]. By using this, we avoid to address the FPT 
problem by Kolmogorov’s equations approach, namely to study the equations associated 
to the two-dimensional process {X{t),Y{t )); many authors (see the references cited above) 
followed this analytical approach to study the distribution and the moments of the FPT 
of integrated BM, and they obtained explicit solutions, in terms of special functions. On 
the contrary, our approach is based on the properties of Brownian motion and continuous 
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martingales and it has the advantage to be quite simple, since the problem is reduced to the 
FPT of a time-changed BM. Actually, for F(0) = ?/ = 0 we present explicit formulae for the 
density and the moments of the FPT of the integrated Gauss-Markov process X, both in the 
one-boundary and two-boundary case; in particular, in the two-boundary case, we are able 
to express the nth order moment of the first-exit time as a series involving only elementary 
functions. 


2 Main Results 

We recall from [1] the following: 

Theorem 2.1 LetY be a Gauss-Markov process of the form fll.ll] .- thenX{t) = x-\- f^^V(s)ds 
is normally distributed with mean x -\- M{t) and variance 7 (p(t)), where M{t) = fom(s)ds, 
j(t) = f^(R(t) — R(s))^ds and R{t) = h 2 {p~^{s))/p'{p~^{s))ds. Moreover, z/ 7 (-|-oo) = 

-|-oo, then there exists a BM B such that X{t) = x + M{t) + B{p(t)), where p(t) = ■j{p(t)). 
Thus, the integrated process X can be represented as a Gauss-Markov process with respect to 
a different BM. 


□ 


Remark 2.2 Notice that, if 7 (-|-cxo) = -|-cxo, though X is represented as a Gauss-Markov 
process for a suitable BM B, X is not Markov with respect to its natural filtration Tt (i.e. the 
cr—field generated by X up to time t). In fact, a Gaussian process X enjoys this property if 
and only if its covariance K{s, t) = cov{X{s),X{t)) satisfies the condition (see e.g. [I5], [30], 
[33]) K{u,t) = u < s < t. Really, if X is e.g. integrated BM with y = 0, x = 0 

(that is, X{t) = J^Bgds), one has K{s,t) = cov(^ B^du, Budu) = ^{3t — s) (see e.g. 
[37] . pg. 654 or [22], pg- 105), and so the above condition does not hold. On the other 
hand, the two-dimensional process ( B^du, B^du) has not the same joint distribution 
as {B{p{s)),B{p{t))), because cov{B{p{s)),B{p{t))) = E[B{p{s)) ■ B{p{t))] = p(s) = s^/3, 
for s < t (see Example 1 below), which is different from K(s, t). However, the process {X, B) 
is Markov, and the marginal distributions of the random vector (X(s),X(t)) are equal to 
the distributions of B{'p{s)) and B{p{t)), respectively; this is enough for the FPT problems 
we aim to investigate. 


Remark 2.3 If 7 (-|-oo) = -l-cxo, and we consider the time average of Y in the interval [0, T], 
i.e. Yt = ? (’ t)y Theorem 12.II we get Y^ = 7 M(T) -|- B{p{T)) , namely, Yt 

is normally distributed with mean {M{T))/T and variance p(T)/T^. In particular, if Y is 
BM, starting from y (that is, m{t) = y, h. 2 (t) = 1, pit) = t), one obtains Yt ^ Xiy, T/3) 
(see Example 1 below and ID)- 


Example 1 (integrated Brownian motion) 

Let be y(t) = y + Bt, then m(t) = y, hi(t) = t, h^it) = 1 and p(t) = t. Moreover, 
R(t) = Jq ds = t and 7 (t) = f^(t — s)^ds = R/3. Thus, p(t) = t^/3, 7 (-|-cxo) = -|-cxo, and so 
there exists a BM B such that W(t) = x -\- yt + 8(^/3) (see [1]). 

□ 
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Example 2 (integrated 0-U process) 

Let Y{t) be the solution of the SDE (Langevin equation): 


dY{t) = -fi{Y{t) - (3)dt + adBt, F(0) = y, 
where p, a > 0 and /S G M. The explicit solution is (see e.g. BD- 

r(i) = /3 + e-'“[!/-/3 + B(p(t)]. 


( 2 . 1 ) 


where B is Brownian motion and p{t) = ^ (e^^* — 1). Thus, F is a Gauss-Markov process 

with m{t) = (3 + e~^^{y — (3), hi{t) = ^ — e~^^) , h 2 {t) = and c(s,t) = hi{s)h 2 {t). 

By calculation, we obtain: 


M{t)= / + e->^\y - P)) ds = pt + 


{y-P) 


m = 




(p )'{.s)ds = 


p 

X _ g-pp"bb 


(1 - e-^^) , 


P 


P ^(s) = TT^ln ( 1 + ) , 


2/i 


cr^ 


7(t) = = 




( 2 . 2 ) 

(2.3) 

(2.4) 

ds 


2a^ 


p? ((j2 + 2pt) fjfi + 2pt/cr2 


P^ Jo \^/l + 2pt/a‘^ -y/rT2/Is7^^ 

2 

^ V^l + 2/it/0-2 - In (l + 2pt/a^) . (2.5) 


mean 


Then, by Theorem 12.11 we get that X{t) = x + f*Y(s)ds is normally distributed with 
x-j-M(t) and variance p(t) = 'y(p(t)). Moreover, as easily seen, limt^+oo 7(t) = +cx), so there 
exists a BM B such that X(t) = x -h M(t) -h B (/o(t)). □ 

Notice that in both Example 1 and 2 it holds p(+cx)) = +cx), so the condition 7 (+oo) = +oo 
is equivalent to p(+cx)) = +oo. 


Example 3 (integrated Brownian bridge) 

For T > 0 and a, /3 G R, let Y(t) be the solution of the SDE: 

dY{t) = dt + dBt, 0<t<T, y(0) = y = a. 

This is a transformed BM with hxed values at each end of the interval [0, T], Y[0) = y = a 
and Y(T) = (3. The explicit solution is (see e.g. [36]): 

/■* 1 

Y{t) =a{l-t/T)+(3t/T+{T-t) / -- dB{s) 

Jo T-s 

= a (1 - t/T) + pt/T + (T - t)B 1 0 < t < T, (2.6) 
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where B is BM. So, for 0 < t < T, y is a Gauss-Markov process with: 


m{t) = a (1 — t/T)+f3t/T, hi{t) = t/T, h 2 {t) = T—f, p{t) 


rp^rj._^y c(sG) = hi{s)h2{t). 


Notice that now p{t) is defined only in [0,T). By calculation, we obtain: 


M{t) = R{t) 


TH{2 + Tt) 

2(1 +Tty ’ 


(2.7) 



T^s y (TH(2 + Tt) T^s(2 + Ts)V ^ 

1+^’ ^ ^ “ Jo V 2(1 + Tf)2 “ 2(1 + Ts)2 ) 


Then, by Theorem 12.11 we get that X(t) = x + f^V(s)ds is normally distributed with mean 
X + M(t) and variance p(t) = 'j(p(t)). As easily seen, limj^r- p(t) = +oo; moreover, by a 
straightforward, but boring calculation, we get that linii^T-- p(t) = 7 i(+oo) = +oo, so there 
exists a BM B such that X(t) = x + M(t) + B (p(t )), f G [0,T]. 

□ 


Example 4 (the integral of a generalized Gauss-Markov process) 

Let us consider the diffusion Y (t) which is the solution of the SDE: 

dY(t) = m'(t)dt + a(Y(t))dBt, Y (0) = y, 

where a{y) > 0 is a smooth deterministic function. In this Example, we denote by p(t) 
the quadratic variation of Y(t), that is, p(t) := {Y)t = Jq a'^(Y(s))ds, and suppose that 
p(-|-cxo) = -l-cxo. By using the Dambis, Dubins-Schwarz Theorem (see e.g. [36]), it follows 
that Y(t) = m(t) + 13(p(t)), t >0 (m(0) = y), where 5 is BM; here, p(t) is increasing, but 
not necessarily deterministic, namely it can be a random function. For this reason, we call Y 
a generalized Gauss-Markov process. Denote by A the “inverse” of the random function p, 
that is, A(t) = inf{s > 0 : p(s) > t}; since p(t) admits derivative and p'(t) = a‘^(Y(t)) > 0, 
also A'(t) exists and A'(t) = we focus on the case when there exist deterministic 

continuous functions a(t), (3(t) (with q;( 0) = /3(0) = 0) and ai(t), /Si(t), such that, for every 
f > 0 : 


a(f), l3(t) are increasing, a(t) < p(t) < j3(t), and ai(t) < A'(t) < /Si(f). 

Since p{t) is not, in general, deterministic, we cannot obtain exactly the distribution of 
Jp Y(s)ds, however we are able to find lower and upper bounds to it. In fact, we have: 

Y Y Y ^ Y fpA) ^ 

/ Y(s)ds= / m{s)ds + / B(p(s))ds = / m(s)ds + / B(v)A'(v)dv. 

Jo Jo Jo Jo Jo 

By using the arguments leading to the proof of Theorem 12.11 (see [1] for more details), we 
conclude that, for fixed t the law of fgY(s)ds, conditional to p(t), is normal with mean 
M(f) = J^m(s)ds and variance 'y(p(t)), which is bounded between ''y(a(t)) and 7 (/ 5 (f)). 
Here, 7 (f) = /J(i?(t) — R(s)yds, where R(t) = f^^A'(s)ds is bounded between ai(s)ds 
and Jq (3i(s)ds. The closer a(t) to (3(t) and ai(t) to (Ji(t), the better the approximation 
above; for instance, if a(y) = 1 -|- ecos^(y), e > 0, we have p(t) = f^(l -|- e cos'^(Y(s))yds 
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and so a(t) = t, (3(t) = (1 + e)^t, ai(t) = 1/(1 + e)^, (3i(t) = 1. The smaller e, the closer 
7 (a(t)) to 

□ 


In the sequel, we suppose that all the assumptions of Theorem 12.11 hold, and 7 (+oo) = 
+cxo; we limit ourselves to consider the special case when m{t) is a constant (that is, m{t) = 
y(0) = y, Vt), thus Y{t) = y + h 2 {t)B{p{t)) and X{t) = x + yt + h 2 {s)B{p{s))ds. Our 

aim is to investigate the FPT problem of X, for one or two boundaries. One approach to the 
FPT problem of X consists in considering the two-dimensional process {X{t),Y{t)) given 
by: 

j X{t) = X + f*Y (s)ds 
|y(t) = + h2{t)B{p{t))dt , 

or, in differential form: 

jdX{t) = Y{t)dt 

\^dY{t) = h'2{t)B{p{t))dt + h2{t)^Jp'{t)dBt , 
and to study the associated Kolmogorov’s equations. 

Many authors (see e.g. |IH], [2l], [2S], [26], |2H]) followed this way in the case of integrated 
BM, namely for Y {t) = y + Bt ■ fact, for r = Tq or r = Ta,b, the law of the couple 
{r{x,y), BT-(^x^y)) was investigated. Let us denote by S the generator of {X,B), that is: 

if one considers, for instance, the one boundary case, then the Laplace transform of 
{Ta{x, y), , dehued for a: < a, ye M, by u{\, v) := E [exp [-\Ta{x, y) - 

(A, z/ > 0), is the solution of the problem with boundary conditions: 

{ ^u{x,y) = \u{x,y), x < a, y eR 
u{a-,y) = y >0 (2.8) 

u{a+,y) = e’^y, y <0 

(see e.g. [26], Lemma 3, or ref. [4], [5], [7], therein). Moreover, for n = 1, 2,... the nth order 
moments Tn{x,y) = E{T^{x,y)) are solutions to the equations = —nT„_i (Tq = 1), 
subjected to certain boundary conditions; however, these boundary value problems are not 
well-posed (see [5T], where some numerical methods to estimate T„ were also considered). 

Notice that, in the case of integrated BM, explicit, rather complicated formulae for the 
joint distribution of [Ta{x,y), BT-^(^x^y)j (and therefore for the density of Ta{x,y)) were found 
in ra. m, H). In order to avoid not convenient formulae, we propose an alternative 
approach, based on the representation of the integrated process X as a Gauss-Markov pro¬ 
cess, with respect to the BM B (see Theorem 12.Ih : this way works very simply, almost in 
the case when y = 0. Thus, in the following, we suppose that Y{t) = y + h 2 (t)B{p{t)) and 
7 (-|-oo) = -|-oo, so the integrated process is of the form X(t) = x + yt + B{p(t)), where 
pit) = 7 (p(f)) and 5 is a suitable BM. Notice however, that the integrated OU process and 
the integrated Brownian bridge belong to this class only ii y = (3 (see fl2.2p ). and a = (3 = y 
(see fl2.7p L respectively. 
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2.1 FPT through one boundary 

Under the previous assumptions, let a be a fixed constant boundary; for a: < a and y G M, 
the FPT of X through a can be written as follows: 

Ta{x, y) = inf{f > Q \ X + yt + B{p{t)) = a}. (2.9) 

Thus, if we set Ta(x,y) = p{Ta{x,y)), we get: 

Ta{x, y) = inf{t >0:5* = h(f)}, (2.10) 

where h{t) = a — x — yp~^(t), and so we reduce to consider the FPT of BM through a curved 
boundary. Since, for x < a and y > 0 the function h(t) is not increasing, we are able to 
conclude that Ta{x,y) is finite with probability one, if y> 0. In fact, as it is well-known, the 
FPT of BM through the constant barrier h(0) = a — x, say f{x), is finite with probability 
one; then, if ?/ > 0, from h(t) < h(0) we get that Ta{x, y) < r(x) and therefore also Ta{x, y) is 
finite with probability one. Finally, if ?/ > 0, we obtain that P{Ta{x,y) < -|-oo) = 1, because 
Ta{x,y) = p~^{Ta{x,y)) < p~^{fa{x)). Note, however, that this argument does not work for 
y <0. 

A more difficult problem is to find the distribution of Ta{x,y), and then that of Ta{x,y). 
However, if h{t) is either convex or concave, then lower and upper bounds to the distribution 
of Ta{x,y) can be obtained by considering a “polygonal approximation” of h(t) by means 
of a piecewise-linear function (see e.g. 0, 0), but in general, it is not possible to find the 
distribution oiTa{x,y) exactly. 


Remark 2.4 Actually, it is possible to find explicitly the density of the FPT of X through 
certain moving boundaries. Indeed, denote by V the family of continuous functions which 
consists of curved boundaries v = v(t), t > 0, n(0) > 0, for which the FPT-density of 
BM through v is explicitly known; this family includes linear boundaries v{t) = at + b (see 
i). quadratic boundaries v{t) = a — bB (see e.g. [7], [20], [38]), square root boundaries 
v{t) = a\/t + b, and v{t) = a a/ (1 -|- bt){l + ct) (see e.g. 0. im. P!). and the so-called 

Daniels boundary v{t) = b — ^ log -|- \l^ (see [13], [E]). For a boundary 

n G V, denote by fv(t\x) the FPT-density of BM starting from x < n(0) through the 
boundary n; if S{t) = v{p(t)) + yt, then one can easily find the density of the FPT of X 
through S, with the condition that x < S'(O) = n(0). In fact, if ts{x, y) = inf{t > 0 : X(f) = 
S'(f)|X(0) = x,y(0) = I/}, one gets Ts{x,y) = inf{t > 0 : x + ty + B{p{t)) = S{t)}; then, 
%{x,y) := p{Ts{x,y)) = inf{f > 0 : x -|- B(t) = v(t)} has density /„ and so the density of 
Ts{x,y) turns out to be 

fsit\x) = f^{p{t)\x)^{t). (2.11) 

For instance, if X is integrated BM (p(f) = t^/3), and we consider the cubic boundary 
S{t) = a + ty + bP (a > 0, 6 < 0), it results S{t) = v{p{t)) + yt, with v{t) = a + 3bt and so, 
for X < a, %{x,y) is the FPT of BM starting from x through the linear boundary a + 3bt. 
Thus, %{x,y) has the inverse Gaussian density fv{t\x) = (see e.g. 0); 

then, the density of Ts{x,y) is obtained by (12.111) . 
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Formula fl2.10p . with |/ = 0, allows to find the density of ra(x, 0) in closed form; in fact, 
Ta{x, 0) is the FPT of BM B through the level a — x > 0, and so its density is: 

/.(iW := 0) < () = (2,12) 

from which the density of ra(x, 0) = p“^(ra(x, 0) follows: 

Ut\x) := |p(r.(x,0) < () = ImWit) = 

If X is integrated BM, we have X{t) = x + B{p{t)), with p{t) = t^/3, so we get (cf. [T7]): 

If X is integrated OU process, the density of Ta{x, 0) can be obtained by inserting in fl2.13p 
the function p(t) deducible from Example 2, but it takes a more complex form. 

Remark 2.5 Formula fl2.13p implies that the nth order moment of the FPT, E(r”(x, 0)), 
is finite if and only if the function is integrable in (0, +cxd). 

Now, let us suppose that there exists a > 0 such that p(t) ~ const ■ t", as t ^ +oo; then, in 
order that E{t^{x, 0)) < oo, it must be a = 2{n + 6), for some 5 > 0. For integrated BM, we 
have a = 3, then for n = 1 the last condition holds with 6 = 1/2, so we obtain the finiteness of 
E{Ta{x, 0)) (notice that the mean FPT of BM through a constant barrier is instead infinite). 
Of course, this is not always the case; in fact, if X is integrated OU process, we have 
pit) const ■ 7 (f) ~ const ■ \n{2pt/a'^), as f —)■ +cxo, and so p{t) = 7 (p(f)) ~ const ■ t, 

as f —)■ +CX 0 , namely a = 1 and the condition above is not satisfied with n = 1; therefore 
E{Taix,0)) = +CX 0 . Not even U((ra(x, 0))^/^) is finite, but U((ra(x, 0))^/^) is so. Notice that 
the moments of any order of the FPT of (non integrated) OU through a constant barrier are 
instead finite. 

As for the second order moment of the FPT of integrated BM, instead, we obtain 
E [(ra(x,0))^] = + 00 , since the equality a = 2(n + 5) with a = 3 and n = 2 is not 
satisfied, for any 5 > 0. 


From fl2.12p we get that the nth order moment of ra(x, 0), if it exists finite, is explicitly 
given by: 

E[iraix,t])Y] = E [(p-^(ra(x,0)))'^] 

"+°° . n-rr . 

(2.15) 






For instance, if X is integrated BM, one has: 


0 + 00 




3^/^(a — x) 


b+cx) 




g-(a-x)V2t^p 


'+CX) 
















By the variable’s change z = 1/t, the integral can be written as: 



1 


{a-xfzl2dz = 


r (1)21/6 r+oo 

(a - x)V3 


(a — xY 


1/6 


. 1 / 6 - 1 ,- 




rQ) 


‘dz 


r (h 2’-''’ 

(a — x)^/3 

where we have used that the last integral equals one, because the integrand is a Gamma 
density. Thus, for integrated BM, we finally obtain: 


E{Ta{x,0)) 




( 2 , 16 ) 


Until now we have supposed that the starting point x < a is given and fixed. We can 
introduce a randomness in the starting point, replacing X(0) = a: with a random variable rj, 
having density g{x) whose support is the interval (—cxo,a); the corresponding FPT problem 
is particularly relevant in contexts such as neuronal modeling, where the reset value of the 
membrane potential is usually unknown (see e.g. Ell). In fact, the quantity of interest 
becomes now the unconditional FPT through the boundary a, that is, inf{t > 0 : X(t) = 
a|y(0) = y}; in particular, if X is integrated BM and y = 0, one gets from fl2.16p that the 
average FPT through the boundary a, over all initial positions rj < a, is: 

T„ = y E{Ta{x,0))g{x)dx = [-] J {a - xY^^g{x)dx. (2.17) 


For instance, suppose that a — rj has Gamma distribution with parameters a, A > 0, namely, 
7] has density 


9{x) 


A" 

fR 


g-A(a-a;) 


(a - 


I(—oo,a) (^) • 


Then, by the change of variable z = a —x one obtains that the above integral is nothing but 
E , where Z is a random variable with the same distribution of a — then, recalling 

the expressions of the moments of the Gamma distribution, one obtains E 
Finally, by inserting this quantity in (I2.17p . it follows that: 


T 


a 


(jfeU.Uillbll. 

y/n F(q:) 


Remark 2.6 For y = U(0) = 0, we have considered the FPT of X through the boundary a 
from “below”, with the condition x = X(0) < a; if one considers the FPT of X through the 
barrier a from “above”, with the condition X(0) > a (namely, inf{t > 0 : X(t) < a|X(0) = 
x, Y (0) = 0}), then in all formulae a — x has to be replaced with x — a. More generally, if 
one considers the first hitting time of X to a (from above or below), a — x must be replaced 
by \a — x\. 
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2.2 FPT in the two-boundary case: first exit time from an interval 

Assume, as always, that 7 (+cxo) = +cxo; for x G (a, b) and |/ G M, the hrst-exit time of X 
from the interval (a, b) is: 

Ta,b{^, y) = inf{t >0 : x + yt + B{p{t)) ^ (a, 6)}. (2.18) 

Set Ta^b{x,y) = p{Ta,b{x,y)), then: 

Ta,bix,y) = ini{t > 0 X + Bt < a-y'^^{t) or x + Bt > b - y'/X^{t)}. (2.19) 

If Ta^b{x,y) is hnite with probability one, also Ta,b{x,y) is so. In the sequel, we will focus on 
the case when y = 0, namely we will consider Ta^b{x,0) = p~^(ja^bix, 0)), where Ta^b{x,0) = 
inf{f > 0 : X + Bt ^ {a,b)}; as it is well-known, Tafi{x,0) is hnite with probability one and 
its moments are solutions of Darling and Siegert’s equations (see [12]). 

First, we will hnd sufficient conditions so that the moments of Tafi{x,0) are hnite; then, 
we will carry on explicit computations of them, in the case of integrated BM. 


Proposition 2.7 If p is convex, then E {Tafi{x,If)) < cxo; moreover, if there exist constants 
c, (5 > 0, such that 0 < p~^(t) < c -1^, then E {Tafi{x, 0))" < oo, for any integer n. 


Proof. If 'p is convex, then is concave, and the hniteness of E{Ta^b{x,0)) follows by 
Jensen’s inequality written for concave functions. Next, denote by f-a,a(t\x) the density of 
the hrst-exit time oi x + B^ from the interval {—a, a), a > 0; we recall from jl2] that the 
Laplace transform of f-a,ait\x), namely, jQ°°e~^^f-a,a{t\x)dt is: 


£ 


cosh(\/^x) 

f-a,a (^ f ) = — , —a < X < a, 6 > 0. 


cosh(\/^aj 

By inverting this Laplace transform, one obtains (see 


( 2 . 20 ) 


f-aAt\x) = ^ (k + 


k=0 


COS 


TTX 

a 


exp 


- k + 


1 \ 2 Oj. 

1 \ X t 


2a2 


( 2 . 21 ) 


The case of an interval (a,b), b > a, is reduced to the previous one; in fact, as easily seen, 
if a = (6 — a)/2 one has: 

fa,b(^M = ’ 

Of course, the density of Ta,b{x, 0) turns out to be fa,b{p{t)\x)fl{t). For the sake of simplicity, 
we take a = —a, b = a, a > 0; then, for x G (—a, a) and an integer n : 


E[{Ta,bix,0))"'] = E [(r_„,a(a:,0))"] = E [(p“^(r_a,„(a:, 0))'"] ='^Ak{x), (2.22) 

k=0 


where 

A,(x) = (k + cos ((i + ^) (r'(i))"*. (2.23) 
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The integral can be written as: 


2^2 

7r2(fc + 1/2)2 


E 


{r\z,))\ 


where Zk is a random variable with exponential density of parameter Xk = (fc + l/2)^7r2/2a^; 
so: 


A(.) = (-1)‘ co= ((i+1) H) ir\z.)r. 


Recalling that E[{Zk)"'^] 
const ■ ; thus: 


by the hypotheses we get E {{p < c^E[{ZkY^] = 


\Ak{x)\ < 


const' 

(fc + l/2)i+2-5’ 


from which it follows that the series J2k ^k{x) is absolutely convergent for every x G (—a, a), 
and therefore E [(r_Q,,a(a;, 0))”] < +cxo. The finiteness of E [(ra,b(a:, 0))"] in the general case 
is easily obtained. 

□ 


Remark 2.8 The condition 0 < p ^(f) < c • is satisfied e.g. for integrated BM, since 

p-l(t) = 31/3^1/3 

(see Example 1), and for integrated OU process, because from the expression 
of p(t) deducible from Example 2, it can be shown that cit < p(t) < C 2 t for suitable ci, C 2 > 0 
which depend on p and a, and therefore < p~^(t) < ^t. 


Now, we carry on explicit computations of E [Tafi{x, 0)] and E [(ra^b(a:, 0))^] , in the case 
of integrated BM. Inserting p{t) = t^/3, {p~^{y) = (Sp)^/^), and n = 1,2 in fl2.22p . fl2.23p . 
after some calculations we obtain: 


E [Ta,b{x, 0 )] 


3V327/3r(|)(6_a)2/3 


TT 


5/3 




/c=0 


{2k + 1)5/3 


COS 


■n{2k + 1) 
b — a 



(2.24) 


A:=0 

Notice that it is arduous enough to express the sums of the Fourier-like series above in 
terms of elementary functions of x G (a, b), and then to obtain the moments of Ta,b{x, 0) in a 
simple closed form; actually, by using the Kolmogorov’s equations approach, in [HT], [32], it 
was obtained a formula for E{Tafi{x,0)) in terms of hypergeometric functions. This kind of 
difficulty does not arise, for instance, in the case of (non-integrated) BM; in fact, by using 
formula fl2.22p with p(t) = t and n = 1, one obtains: 


(2fc + l)4 


cos 


7r(2/c + 1) / a + b 

—j - I X - 

b — a 


(2.25) 


E [r-aA^)] = 


32a2 






k=0 


(2fc + 1)3 


COS 


71 


(2k + l)-x 


on the other hand, the well-known formula for the mean first-exit time of BM from the 
interval {—a, a), provides that the sum of the series must be — x‘^. 
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Figure 1: Plots of the mean exit time, 0)), of integrated BM from the interval 

(—1,1) (lower curve), and of the function z{x) = 1.35- (1 — (upper curve), as functions 

of x e (—1,1). 


However, fl2.24p and fl2.25p turn out to be very convenient to estimate the first two 
moments of Ta,b{x, 0) for integrated BM; in fact the two series converge fast enough, so to 
obtain “good” estimates of the moments, it suffices to consider a few terms of them. As for 
E [Ta,b{,x, 0)], it appears to be htted very well by the square root of a quadratic function. 
In the Figure 1, for integrated BM, we compare the graphs of E{Ta,b{,x,0)), calculated by 
replacing the series in (I2.24p with a finite summation over the first 20 addends, and that of 
C ■ [{b — x){x — a)]^/^, as functions of x G (a, 6), for a = —1, = 1, and C = 1.35; the two 

curves appear to be almost undistinguishable. 

We have also calculated the second order moment of the first-exit time of integrated 
BM, by summing the hrst 20 addends of the series in fl2.25p . In the Figure 2, we plot 
E [(ra^b(x, 0))^] , E"^ [ra^b(x, 0)] and the variance 

Var [ra,6(x, 0)] = E [(ra,b( x,0))^] — {E [Ta,b{x,0)])‘^, as a function of x G (—1,1), for a = 
—1, 6 = 1; as we see, the maximum of Var [ra^b(x,0)] is about 10% times the maximum of 
F;(r_i,i(x, 0)). 

As in the one boundary case, if we introduce a randomness in the starting point, replacing 
A(0) = X G (a, b) with a random variable r], having density g{x) whose support is the interval 
(a, b), we can consider the average exit time over all initial positions r; G {a,b). li y = 0, this 
quantity is: 

Ta,b= f E{Ta,b{x,0))g{x)dx. 

J a 

In the case of integrated BM, Ta,b can be calculated by using the expression of E{Ta,b{x, 0)) 
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Figure 2: From top to bottom; plot of the second moment (first curve), the square of the 
hrst moment (second curve), and the variance of the first-exit time r_i^i(x,0) (third curve) 
of integrated BM from the interval (—1,1), as functions of x G (—1.1). 


given by fl2.24p . We obtain: 

3l/327/3rM)(^_^)2/3 


Ta,b — 


yj-5/3 




A:=0 


{2k + 1)5/3 


7i{2k + 1) 
b — a 


a + b 


cos ^ I X -::— I g{x)dx 

(2.26) 

(it has been possible to exchange the integral of the sum with the sum of the integrals, 
thanks to the dominated convergence theorem); the integral in fl2.26p equals E{Uk), where 


Uk = cos 


7r(2/c+l) 

b—a 




< 1. Therefore: 


^ 3‘/W(|)(6 - aps 

-■ «.!' --=575- 


TT 


k=0 


(2fc + 1)5/3 


E{Uk). 


(2.27) 


In the special case when g is the uniform density in the interval {a,b), we get by calcu¬ 
lation: 


■' «.»-=575- 


yj-5/3 


{2k + 1)5/3 


COS 


fc =0 

3l/32l0/3r(4)^^_^)2/3 

);t73 


7r(2/c -|- 1) 


a + b 


X 


dx 


E 

fe =0 


{2k + 1)8/3 • 


(2.28) 


Thus, Tafi = const - {b — a)'^^^. This conhrms the result by Masoliver and Porra (see [3T], [32]). 
obtained by the Kolmogorov’s equations approach in the case of integrated BM, with y = 0 
and uniform distribution of the X— starting point, according to which, the dependence of 
Ta,b on the size L = {b — a) of the interval, is L^/3_ 

As far as integrated OU process is concerned, the moments of Tafi{x, 0) can be found again 
by formula fl2.22p . where p(t) can be deduced from Example 2; however, it is not possible 
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Figure 3: Plot of numerical evaluation of the mean exit time, E (r_i 4 (x, 0)), of integrated 
OU with = 1 / = 0, from the interval (—1,1), as a function of x G (—1,1), for a = 1 
and several values of fx. From top to bottom, with respect to the peak of the curve: /i = 
2; 1.8; 1.6; 1.4; 1.2; 1. 

to calculate explicitly the integral which appears in the expression of Ak{x), so it has to be 
numerically computed. Since the integrand function decreases exponentially fast, it suffices 
to calculate the integral over the interval (0,10), to obtain precise enough estimates. In the 
Figure 3 we have plotted, for comparison, the numerical evaluation of the mean exit time of 
integrated OU process with y = f3 = 0, from the interval (—1,1), as a function of x G (—1,1), 
for a = 1 and several values of p; in the Figure 4 we we have plotted the numerical evaluation 
of E (x,0))^] , [r_i_i(x, 0)] and the variance Var [r_i^i(x,0)] of the hrst exit time 

of integrated OU process, for a = 1 and /i = 1. As we see, the maximum of Var [r_i_i(x, 0)] 
is about 5% times the maximum of E (r_i,i(x, 0)). 

Finally, we mention the exit probabilities of the integrated Gauss-Markov process X 
through the ends of the interval (a, b), namely: 

TTa{x,y) = P{Taix,y) < n{x,y)) = P {X{Ta,bix,y)) = a ), 

and 

TTb{x,y) = P{Tb{x,y) < Taix,y)) = P {X{Ta,bix,y)) = b ). 

Recalling the well-known formulae for exit probabilities of BM, we get, for ?/ = 0 and 
X G (a, b) : 

7ra(x, 0) = P (^x -h B{Ta,b{x, 0)) = oj 7rfe(x, 0) = P ^x -h B{Ta,b{x, 0)) = 6j 

Notice that, in the case of integrated BM, several probability laws related to the couple 
[xa^b, J were evaluated in [21] (in particular, explicit formulae for 7ra(x,0) and 7rfe(x, 0) 
were obtained), but they are written in terms of special functions. 
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Figure 4: From top to bottom: plot of the second moment (first curve), the square of the 
hrst moment (second curve), and the variance of the hrst-exit time r_i^i(x,0) (third curve) 
of integrated OU with ?/ = /3 = 0, from the interval (—1,1), as a function of x e (—1,1), for 

(T = 1, /i = 1. 
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